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ABSTRACT 


__ A plasma resonator experiment gives several resonances, whereas theory predicts only one 
for a homogeneous zero-temperature plasma cylinder of circular cross section. A theory is 


presented that predicts: that also an inhomogeneous plasma cylinder generally has only one 


dipole resonance; that a deformation of the cross section or introduction of a longitudinal 
“magnetic field splits the resonance into two; and that an elementary model of a plasma of 


non-zero temperature gives an infinite series of resonances. This theory agrees qualitatively 


/with the experimental results. 


The surface on which the dielectric constant vanishes in an inhomogeneous plasma consti- 
tutes a mathematical singularity. A suitable treatment is presented that resolves the paradox 
of losses in a plasma where, by assumption, there are no collisions. 


1. Introduction 


As early as 1931 Tonks [1] put a low pressure mercury discharge tube of circular 
cross section between two condenser plates loading a Lecher line. By varying the 
current density in the tube he found one strong resonance and a weaker one. Several 
years later Shapiro [2] reported measurements of a similar character. He, however, 
found only one resonance. These two experiments were done at a fairly low frequency, 
a few hundred Mc/s. 

In connection with his studies of the scattering of electromagnetic waves from 
meteor trails, Herlofson [3] realized the value of laboratory experiments. In 1951 
Romell [4] thus published results of some measurements made at this laboratory 
on the scattering of electromagnetic waves from a low pressure mercury discharge 
tube. The electric vector was perpendicular to the axis of the tube. The frequency he 
used was 109 c/s, i.e. a frequency a few times greater than that used in the two experi- 
ments just mentioned. Romell reported three resonances with decreasing amplitude 
for decreasing plasma density. 

Then, in 1957, Dattner [5] published the results of measurements on the energy 
transmitted through a wave guide containing a discharge tube. The electric field in 
the wave guide was perpendicular to the axis of the tube. Dattner observed minima 
in transmitted energy at approximately the same densities as those at which Romell 
had observed maxima in reflected energy. The frequency was five times that used 
by Romell, and the diameter of the discharge tube was decreased in the same ratio. 
This should make the effect of the tube the same as in Romell’s case. 
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Romell measured the space field distribution round the tube and found that the 
scattered field was a dipole field. Several years later Boley [6] repeated the same 
measurement and came to the same conclusion for the two major resonances. Even 
if the measurements are conclusive only for the major resonances, it seems quite 
plausible that all the observed resonances have dipole space distribution. 

Tonks reported measurements with a discharge tube of elliptic cross section. From 
his curves one finds that there is one resonance corresponding to each of the axes 
of the cross section (Tonk’s conclusion is somewhat different). In Dattner’s experiment 
we also find that each single resonance for circular cross section splits into two reso- 
nances for elliptic cross section—corresponding to the two axes of the ellipse. 


he 
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In simple resonance systems we generally find two kinds of energy which we may 
call potential energy and kinetic energy. A plasma cylinder is no exception. In the 
simplest case, which we shall deal with here, the plasma density is not a function of 
the coordinate along the axis of the cylinder. Let us begin with a cylinder where the 
space charge density is zero everywhere and all particles are at rest. If we then perturb 
the position of the electrons, an electric field is created and this field tends to decrease 
the perturbation. The electrons are set into motion. As usual, when the perturbation 
in position has become zero, the particles have a non-zero velocity and swing out 
in the direction opposite to the original perturbation. Again we find that a field is 
created that tends to diminish the perturbation. The system is oscillatory. The two 
kinds of energy we find here are electric field energy, which has to be regarded as 
potential energy, and kinetic energy of the electrons. The positive ions and possible 
negative ions present are so heavy that, near the resonance frequency, their motion 
can be neglected and consequently also their contribution to the kinetic energy. Since 
the energy of the system is periodically transformed from kinetic energy into potential 
energy and back again, the condition for resonance is that the total energy is inde- 
pendent of time. Since the potential energy is not restricted to the interior of the cyl- 
inder, the resonance frequency generally depends not only on the conditions inside 
the plasma region, but also on the shape of the plasma region and the conditions out- 
side. 

For every particle in the medium we have the equation of motion 


tomv=eH, 


where v is the velocity of the particle, e and m its charge and mass, and E the electric 
field. If the perturbation is periodic with the angular frequency w, we then easily 
find that the kinetic energy per unit volume is 
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where n is the electron density. If we denote the electric field energy per unit volume 
by W, we find that the dielectric constant, as introduced by Bergmann and Diiring 
1929 [7], can also be written 
a) cep 
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We stress that 7 and W are time averages of the energy densities and not of the total 
values for the whole system. 
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generally there are regions which do not contain plasma, i.e. regions where 
netic energy is zero, we must also have regions where the kinetic energy exceeds ° 
potential energy, as equality between total kinetic energy and total potential 
‘gy - necessary for resonance. This means that there must be regions where ¢ 
egative. re a 
Actually we find that, for resonance, the integral of 
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Daxtsnitet over the whole space must vanish. (At this point it is perhaps necessary 
_ to stress that, in systems where the fields are independent of time, the above quantity 
represents the total potential energy of the system. For time-dependent fields this 
_ is approximately correct at frequencies far below the resonant frequencies of the 
_ system; but near resonance the kinetic energy of the system must be taken into 
_ account and then the integral is no longer equal to the energy of the system but to 
_ the difference between the potential and kinetic energy.) 

In the case of a cylinder of uniform density the electric field does not vanish outside 
the cylinder. Since the kinetic and potential energy must be equal at resonance we 
_ find that ¢ for the plasma region must be negative. In this special case the electric 
_ field energy is equally divided between the inner and the outer region if we assume 

static field distribution; we find that resonance occurs when ¢ = — 1, a condition 

given by Tonks [1]. 
Up to now we have assumed that the only two kinds of energy of interest are the 
kinetic energy of the electrons in the plasma region and the electrostatic field energy. 
_ This means that we have neglected the magnetic field energy, which is acceptable if 
the extension of the plasma region is small enough. We have also neglected the thermal 
motion in the plasma. With these limitations we have found that we cannot expect 
any resonance effects if the dielectric constant is positive everywhere in space. 

As ¢ also has the ordinary meaning, the field distribution can be determined in the 
usual way. For any volume consisting of a finite number of homogeneous regions we 
will always find finite values for the field quantities—also in the case where ¢ = 0 in 
some region or regions. In particular, we may mention that any region enclosed by a 
region with ¢=0 is completely shielded and that the potential is constant on the 
surface of the shielded region. This means that the normal value of the electric field 

' strength at the inner surface may increase indefinitely with decreasing thickness of 
the enclosing shell. This discussion indicates that the surface ¢ = 0 in an inhomo- 
geneous medium may be a singular surface, on which the electric field strength is in- 
finite. Since infinities do not exist in nature there must be some unphysical assump- 
tion in our reasoning. Let us therefore for a moment study the concept of the dielectric 
constant of a plasma. 

True continuous space charge does not exist in reality since charge is quantized 
into particles. For a plasma we thus ought to introduce the charges as singularities 
in vacuum. Because of the enormous number of singularities it is usually convenient 
to average the field quantities over regions that are large compared with the inter- 
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particle distances and in that way be able to treat the plasma as a continuous medium. 
This point of view led to the concept of the dielectric constant of a plasma. With 
reference to Lorentz’s work [9] on the dielectric constant of solid and liquid media, 
Hartree [10] criticized the derivation of the plasma dielectric constant on the ground 
‘that the Lorentz polarization term was not included [9]. In two important papers 
Darwin has then proved that the original definition is the correct one in the case of 
a plasma region [11, 12]. 

The continuum treatment should be valid when the relative variation of the field 
quantities within an interparticle distance is insignificant. In regions where ¢ is differ- 
ent from zero this condition is generally fulfilled but in regions where e changes sign 
by passing through zero a careful study is necessary. To exemplify, we can look at a 
plane condenser where the electric displacement, D, has the same value everywhere 
between the condenser plates. If ¢ varies linearly with the coordinate perpendicular 
to the plates and passes through zero at some plane, the electric field, H, must become 
infinite as the inverse of the distance to this plane. Then there must be a region 
near that plane in which the relative variation of the electric field over an interparticle 
distance is of the order of unity or larger. This clearly shows that a real medium cannot 
be treated as a continuous medium in the vicinity of a region where ¢ is zero (cp. [13], 
[14]). 

_In the case of an inhomogeneous plasma it is difficult to obtain the solution in a 
closed form using tabulated functions and to find the solutions by numerical methods 
would be very laborious. Therefore we may ask whether it is possible to approximate 
the continuous distribution by consecutive homogeneous regions. To gain some insight 
into this problem let us for a moment return to the plane condenser and let us assume 
that the dielectric medium consists of a number of homogeneous plasma regions of 
different density. The electric displacement, D, must then have the same value in all 
regions. When the frequency is such that ¢ is almost zero in one of these regions, the 
electric field is much larger in that region than in any other of the layers. For a 
given value of D, the difference between the potential and the kinetic energy in our 
dielectric medium tends to plus infinity when ¢ tends to zero from one side and to 
minus infinity when we approach from the other side. Thus between two consecutive 
singular values there must be some intermediate frequency for which this difference 
is zero and we should expect as many resonances as there are layers. This simple 
example seems to indicate that the density in an inhomogeneous region cannot be 
approximated by a step function. 

Up to this point we have neglected the thermal motion. As a first attempt to inves- 
tigate the effect of a non-zero temperature it may seem reasonable to introduce the 
electron pressure in the macroscopic equation of motion. Doing so one finds that, 
for a given exciting field, the solution is the same as for zero temperature. However, 
one also gets a spectrum of discrete modes representing oscillations inside the plasma 
region, which are not coupled to the field in the space outside. If w is given, there are 
only a finite number of density value at which these oscillations exist. All these values 
fall in the region e > 0. 

To treat the effect of the thermal motion by macroscopically introducing the elec- 
tron pressure is surely a crude method to find the effect of a non-zero temperature. 
A more satisfactory way would be to study the motion of every separate particle 
and to sum up the effects of the individual particles. This method, though compli- 


a is surely the only method that can yield the detailed behaviour of a hot 
plasma. 
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The theoretical discussion of the electromagnetic behaviour of a plasma began 
_ shortly after the electron was discovered. First one investigated the propagation of 
_ electromagnetic waves. Later one became interested in oscillations in plasma regions 

much smaller than the vacuum wave length. J. J. Thomson (1928) [8] may have been 
; the first to publish an investigation in this field, which was followed by one by Tonks 
_ and Langmuir [15]. Three years later, Tonks [1] discussed the influence of the shape 
__ of the plasma region. 

The first reports were restricted to homogeneous plasma region. In 1931 Tonks 

_ also considered an inhomogeneous region in the plane case. He gave only the reactive 

_ component of the impedance, however, and thus avoided the difficulty due to the 

_ singular surface at ¢=0. In the cylindrical case there are two singularities, one at 

_ the surface, e = 0, and, as is usual in cylindrical problems, one at the origin. Generally, 
it is impossible to find a solution that is regular at both singularities, and one has to 
find an artifice which gives an acceptable solution. One then ought to take the solution 
that is regular at the origin. The other method of removing the other singularity has 
been to introduce losses in the plasma (Herlofsson, 1951 [3], Kaiser and Closs, 1952 

[{16]). Having connected the solution on the two sides of the singularity, one lets 

_ the collision frequency tend to zero and finds that, even in the limit of zero collision 
frequency, energy is absorbed. In the case of reflections of plane electromagnetic 
waves the same phenomenon is found (Budden, 1955 [17]). Budden states that the 
group velocity decreases with the distance to the singular surface, to become zero 
there. Hence energy will pile up at this surface. In this way Budden tries to explain 
the difference between the incident and the reflected energy. 

In the case of an inhomogeneous plasma it is difficult to obtain a solution in closed 
form using tabulated functions. Therefore Makinson and Slade [18] tried to approxi- 
mate the continuous distribution by a finite number of homogeneous regions. In the 
case of a homogeneous cylinder there is only one dipole resonance. Makinson and Slade 
found that, for an increasing number of layers, the number of resonances increase, 
but not as rapidly as the number of layers. In this way they tried to explain the series 
of dipole resonances observed by Romell. 

W. O. Schumann [19] attempted a different method to explain the many observed 
dipole resonances. He introduced a plasma region as dielectric medium for a short 
length of a two-wire line. The plasma was inhomogeneous with a discontinuous drop 
in density at the boundary. He assumed a non-zero collision frequency. In this way 
he found, however, only two maxima in the energy passing the loaded part of the line. 

At the Linde conference in 1959, Gould [20] presented a paper where he had intro- 
duced a non-zero electron pressure in the macroscopic equation of motion. By changing 
the electron density from zero to infinity, one finds a finite number of resonances. 
The number of resonances decreases with increasing temperature. Consequently, the 
relative positions of the resonances depend on the temperature. 

Hershberger’[21] also assumes a non-zero temperature. He finds an infinite number 
of maxima and minima in the reflected wave which are due to some kind of inter- 
ference. In contrast with the result obtained by Gould, the relative positions of his 
maxima are independent of the temperature and the tube diameter. 


The present investigation originated with the wish to explain the many dipole 
resonances observed by Romell. Obviously a homogeneous cylinder of circular cross 
section is not a satisfactory model. Merely assuming a different shape of the cross 
section is not enough. We find only a doubling of the number of resonances for an 


167 


i ae 
7 


et ee ve — 2 ee: 


E, AstROM, Electromagnetic behaviour of a plasma 


elliptic cross section and the same result is to be expected for any shape. It thus 
seems that the experimental results cannot be explained on the model of a cold, 
homogeneous region. 

A cold inhomogeneous plasma region exhibits a singularity at the surface e =0, 
which requires attention. According to the usual method one assumes that the pertur- 
bations are so small as to make non-linear effects negligible, that there is no thermal 
motion, and that there are no collisions—and yet there are losses. That fact ought 
to make us suspicious. Actually we find that near and at the singular surface the plasma 
does not behave as a continuous medium, Taking this into account we again investigate 
the scattering from an inhomogeneous plasma cylinder. We now find only one reso- 
nance of dipole character. Further there is no absorption of energy, and the energy 
radiated at resonance is the same as for a homogeneous cylinder of the same di- 
ameter. (G. Boyd at California Institute of Technology has carried out a similar investi- 
gation, but he treated the singularity in the common way.) Although we have investi- 
gated only a special case our result seems to indicate that an inhomogeneity in the 
plasma cannot increase the number of resonances over the number we find for a 
homogeneous plasma region. 

We also investigate the possibility of approximating the continuous distribution 
by a series of homogeneous layers. From a study of a plane condenser with a dielectric 
material of such a character we find that the continuous density distribution cannot 
be approximated by a discontinuous one. 

Since the attempts to explain the dipole resonances on the model of a cold plasma 
have failed, we try to take the effect of thermal motion into account. In order to 
simplify the problem as much as possible, we discuss the plane case with only one 
value on thermal velocity in the unperturbed case. We then obtain an infinite series 
of resonances. At the present stage of development, it seems hopeful that we shall 
be able to explain the experimental results by introducing a non-zero temperature, 
even if this goal has not yet been reached. 


2. Plane geometry 


2.1 Zero temperature 
(a) Introduction 


We begin with the simplest possible case and step by step introduce new features 
and discuss the consequences of each step. In order not to stray too far from the 
main subject we have found it necessary to exclude many interesting questions. For 
instance we have only touched on the question of the stability of a plasma, which 
is such a vast subject that it would require a separate investigation. 

It is natural to begin with the one dimensional (plane) problem. When we then 
introduce the inhomogeneous plasma we come across the singularity at e =0 and we 
have to resolve the difficulty thereby introduced. Our discussion shows that the con- 
tinuity of the density distribution is essential: an attempt to approximate a continuous 
density distribution by a discontinuous one results in a solution that is not a valied 
approximation to the solution of the continuous case. 

For those familiar with circuit theory the introduction of an equivalent circuit will 
prove helpful. Therefore, we have in this paper used that possibility, despite the fact 
that in one case we obtain a capacitance of negative sign. 
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The plasma can also be thought of as the dielectric medium of a two wire (or two 
pee line. The field of research is too large and we have only hinted at this possi- 

} ility. & , 

a 


_ Let us introduce two parallel plane plates of infinite extension. The distance be- 
_ tween them may be d and the potential difference U. The upper one may be defined 
_ to be positive. The medium in between is constituted of free charged particles of some 
_ tew different species. For each species the density may be a function of the coordinate 
_ perpendicular to the plane of the plates but independent of the other cartesian coor- 
_ dinates. We assume that, in the unperturbed state, the particle velocity in the direc- 
_ tion perpendicular to the plates shall vanish. Since we shall not assume that we gener- 
_ ally have space charge neutrality, we instead assume a homogeneous magnetic field 
_ in the direction of the plane of the plates to ensure that an equilibrium is possible in 

the unperturbed state. 

A continuous distribution of space charge never occurs in nature. Space charge is 

_ formed by a more or less regular pattern of discrete charges. For our calculations we 
can, however, with few exceptions, think of the charge as smeared out to form a 

continuous distribution. As this will simplify the mathematical treatment considerably, 

we shall in the following assume that we have a continuous distribution of space 

charge. Since the value of each discrete charge is small and the number of particles 

per unit volume large, the error introduced is generally very small, but we must 
ascertain in every specific example that the smearing out is a suitable approxima- 

tion. 

To simplify the treatment we shall also assume that the distance between the plates 

is much smaller than the vacuum wavelength of the perturbing field. 

In order to be able to calculate the perturbation of the space charge we need to 
know the motion of the individual particles. We then must know the electric field 
along the path of the particle. This field will be different for the different kinds of 
particles present at any given position. 

The electric field that a particle feels may conveniently be divided into three parts. 
The source of the first is the applied voltage U between the plates. This field strength 
equals —U/d. Consequently this field has the same value everywhere between the 

lates. 
P The source of the second part is the displacement of the charges from their equili- 
brium position. As we will find the field created by this source has, in the plane case, 
the same value everywhere between the plates and the field created is independent 
of the position of its source. 

The third source is the relative motion of the different kinds of particles. We here 
count those particles that pass the level at which the test particle, which may be any 
of the particles’ present, is situated. Consequently this part, and this part alone, does 
depend on the condition of the medium at the point where the particle under study 
is situated, and it exists only if there are more than one kind of particle present. 

For the study of the second part we begin by looking at the space charge due to 
one species of particles only. We may then add the similar contributions from other 
kinds of particles. The space charge density we call 9. Let us take a layer of thickness 
dy at the distance y above the lower plate, as shown in Fig. 1. If the voltage between 
the plates is zero we find that the voltage at the position of the charge, due to the 


charge ody, is 
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Fig. 1. An infinitesimal plasma layer between parallel metal plates. 


a ¥ld=y) ody - (1) 
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Consequently the corresponding electric field is 
(1) above the charged region 


nay 27, @) 
d & 
(2) below the charged region 
B= wed eoy (3) 
Os aes 


In the present investigation we are not specifically interested in this field but more 
in the change created by perturbations of the position of space charge. Assume the 
perturbed position to be y + 7(y). The perturbation of the field due to the perturbation 
7 in the position we may call 6#. 


=1 oy 4 
OH 7 , (4) 


We thus find that the field 6H has the same value everywhere between the plates. 
We also find that the change is independent of the position of its source. 

We now come to the third part. We assume a background charge density o. If our 
test particle moves upwards the distance Ay relative to the particles constituting 
the background, then the test particle has passed the charge oAy. To begin with 
this charge layer was above the test particle and the field created was that given by 
equation (3). After the perturbation is introduced the same charge is below the test 
particle and the field felt is that given by equation (2). The change in field due to the 
perturbation of the position of the test particle is thus the difference between the 
values given by the equations (3) and (2). Let us call the perturbation of the position 
of the test particle ,. In these-equations we then have to introduce AY =). — ts 
where uu, is the perturbation of the background charges. The field due to this source 
is thus 


yee (Nu —p)- (5) 


0 
To summarize our result, let us introduce n kinds of particles with densities 
Q,(y) and with perturbed positions y+7,(y). The field E,(y) felt by a test 
particle of the u™ species is then 
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_ if the y coordinate is perpendicular to the plates and if the magnetic field is directed. 
along the z-axis. The first of these equations we can integrate once and obtain 


es eB 
$ | an, (9) 


4 Introducing this into the second equation we find 


} d*7 eB\? eH 
; —+ | — =—, 10 
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When we have many kinds of particles present, we get a considerable simp- 
_ lification if we, instead of 7, use 


vaya hee (11) 


r 
7 
Cu 


. If we now introduce the equations (6) and (11) into the equation of motion (10) 
and write it specifically for the w™ kind of particle, we obtain 


: 2 2 n aig» Y Dunls UF 
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dt? My y=1 M, Eq My Eo M, Eq 
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Thus Q is the net charge density. It vanishes for a neutral medium. 

For the special case that the exciting voltage is harmonic in time, we may introduce 
the admittance Y per unit area of the plates, since we have assumed the distance 
between the plates to be much smaller than the vacuum wavelength for the exciting 
frequency. The charge density on the upper plate is equal to the electric displacement 
at the surface, with opposite sign. The current is the time derivative of that charge 
density. Hence the admittance is 


ty acre ( S | vs ods) Pres (13) 


y=1 


where the first term represents the vacuum capacitance between the plates. The mo- 
tion of the space charge thus represents an admittance in parallel with this vacuum 
capacity. 

If we see the structure as a parallel line for propagation and restrict the discussion 
to those cases where the wavelength is much larger than the distance between the 
plates, we find, as usual, that the phase velocity v, is given by 


“2 = ( tom \% 
c eh} * 


where c is the velocity of light in vacuum. 


(b) One-species medium (electrons) 


In this case we can drop the index. Equation (12) is then simplified to 


d?r (/eB\2 r €0 U 
Tat (S] r= [522 as. (14) 


The right member is independent of the space coordinate and a function of time only. 
Since the coefficients in the left member are also independent of the Space coordinate, 
so is r, Consequently r can be put outside the integral sign. The integral that remains 
is then the mean value of e o/m ¢, between the plates. We wish to stress that also those 
regions that have no space charge must be included when we calculate the mean value. 
If the mean is indicated by a bar we find 


d?r [/eB\? ee U 
dt? | = M Eo = ae (15) 


Let us first assume that the space charge is concentrated to a finite interval in the y 
direction and that the plates are at infinite distance from that region. We then find 
that the charge does not enter our equation at all. The perturbation is the same 
for all particles in the charged region, which means that the Space charge moves as 


3 
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a rigid body. This rigid body can oscillate and the frequency of oscillation is the gyro - 


frequency of the charged particles constituting the space charge. We have thus found 
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I ‘ig. 2. The admittance of uncompensated electron space charge 
mn a magnetic field for the plane case. 


Ley 1 ep [fes\P ent? ae 
fee ee (:2) eh ane 
ie Eo Ed me |\m M Eo M Ey 
Z 


that the system has an equilibrium position, that this equilibrium is stable for pertur- 
_bations of the kind we have discussed, and that this is the case whatever the charge 
distribution within the charged region. 
Now if the distance between the plates is finite, the space charge does influence the 
phenomenon, though only through the mean value of the charge between the plates. 
_ This means that, for a given value of the distance between the plates, our phenomenon 
_ does depend on the total charge between the plates but is independent of the distribu- 
tion of the charge. Also in this case, the perturbation has the same value for all 
particles within the charged region, which means that the charged region moves as 
-a rigid body. Also here we have a natural frequency of oscillation. Its value is smaller 
than in our previous case by an amount given by the mean density of the space 
charge. 
_ When the total charge between the plates exceeds a certain critical value, the solu- 
tion of the homogeneous equation is exponential with real values of the exponent. 
One of these value is positive. Consequently our equilibrium is unstable. Also in 
_ this case, we can find oscillatory solutions of the inhomogeneous equation. Since the 
- equilibrium is unstable the equilibrium cannot be realized in an experiment and thus 
‘the solution of the inhomogeneous equation when the total charge between the 
plates exceeds the critical value is of theoretical interest only. 


The admittance for the case when only electrons are present. The admittance is 
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If we represent the admittance by passive elements, the result will be that given in 
Fig. 2. C is the capacitance when no charges are present. The space charge acts as 
a series resonance circuit where the capacitance is positive for a space charge less 
than the critical value but has a negative value when the space charge exceeds the 
critical value. Because of the dependence of the corresponding admittance on the 

frequency it must, despite the sign, be regarded as a capacitance. Here we must 
also remember that the case of negative capacitance cannot be realized in an expert- 
ment because the plasma is not stable in that case. 


(c) Introduction of a stationary positive background 


Also in this case we would like to have r independent of the y coordinate. We are 
again free to'choose any distribution of the electron space charge density. The positive 
space charge, on the other hand, must be homogeneous in the whole region where 
electrons are present. This different influence of the positive and the negative space 
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charge density is seen clearly in equation (12), where the negative space cha ge 
appears only under the integral in the right member. This means that only the mean 
value of the negative space charge over the whole region between the plates is of 
interest. On the other hand, the local positive space charge density enters in the left 
member. Therefore the negative space charge moves as a rigid body only if the positive 
space charge is homogeneous. 

With the assumption we have made, equation (12) has the form 


r+ {(2)'-ce: <2 ras. a7) 


dt? m _MEy MEo 


Here ¢ and m refer to the electron, @_ is the mean value of space charge density of the 
electrons and g, is the local value of the space charge density of the psoitive ions. 
In the third term e is the charge of an electron and 9, the charge density of the posi- 
tive ions. Since these two quantities have opposite sign, the third term in equation 
(17) is positive. This term corresponds to, and is generally known as, the plasma fre- 
quency. We thus see that the plasma frequency is given by the charge of the moving 
particle and the density of the background space charge. The density of the moving 
species does not enter in this connection. The last term in the left member refers to 


‘the mean value of the electron space charge. It has opposite sign to the second and 


third terms. 

If the charges are confined to a finite interval in the y direction but the plates are 
at infinite distance, the negative space charge density does not enter into equation 
(17). Our system then has a natural frequency which is the root of the sum of the square 
of the gyro frequency of the electrons and the so-called plasma frequency defined in 
the way just discussed. In this case our system is always stable. 

If the distance between the plates is finite, the mean value of the electron space 
charge enters, decreasing the natural frequency of the system. When it exceeds a 
critical value, the equilibrium of our system is unstable as in our previous example. 

In the special case when we have neutrality at every point between the two plates, 
we can have equilibrium also in the case when the magnetic field vanishes. In this 
case the natural angular frequency of our system is 


(: od (18) 


Eo 


If the space charge fills only a very small part of the region between the plates, it is 


equal to the so-called plasma frequency but if the space charge fills the whole region 
it is zero. 


(d) Same case with inhomogeneous background charge density 


In the last section we had two species of particles, mobile electrons and stationary 
ions. We found no difficulty whatever. The space charge distribution of the electrons 
was a function of the y coordinate, but we had to assume that the background charge 
was homogeneous. In this section we shall drop this limitation and allow any distri- 
bution of the stationary background charge density with the y coordinate. 
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. In the equation of motion for the electrons (12), we find that the right member 
is independent of the space coordinate. Let us denote it by S. If we then restrict our 
Bee cussion to the case where r and U are proportional to e'® we find 

Z 


Te 2 -1 
™m M Ey 
i 


where e and m refer to the electron but 9, is the positive background space charge 
density. The latter quantity is a function of the space coordinate and hence also r. 
If we now introduce our expression for r in the right member of equation (12) and 
remember that S does not depend on the space coordinate, this right member be- 
comes 


=. 


2 -1 
Se= 2 with o= 4 | cea (=) eee wt, ds. (20) 
d d M Eq m M Eo 


But the right member, i.e. this expression, is equal to S. Thus we find 


S= 5 (l—o)"*. (21) 


Solutions of the homogeneous equation. In order to get a solution when U vanishes, 
_we must require that the denominator of expression (21) vanishes. Let us assume 
that w? is an arbitrary complex quantity 


w=pt+igq. (22) 


Here p and q are arbitrary real quantities. The imaginary part of the denominator of 


(21) is then 
x ; ate 
~if { 2 |} (<2)'- 222-»} +q ag (23) 
d J méy m ME 


The integrand retains one and the same sign throughout the interval of integration. 
Expression (23) therefore can vanish only if g = 0. Hence the denominator of expres- 
sion (21) can vanish only if w? is a real quantity. This means that values of q@ permit- 
ting solutions of the homogeneous equation must be either real or purely imaginary. 

The case of a parabolic distribution of the background space charge density. As an 
example we may study the case when both the positive and negative space charge 
densities have a parabolic distribution in the y coordinate. We assume the density 
to be zero on the boundary of the space charge region. 


2 
tpt a2 (1 dee a when —a<y<a _ and zero elsewhere, | 
ME of 
(24) 
2 
pera Sate ( = ) when —a<y<a and zero elsewhere. 
ME oe 
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In the second case our expression for 7 becomes infinite for one value of y. Let us 
assume that the plasma density and the space derivative are continuous functions 
of the space coordinate in a region which includes this point. Then, if y =Y, is the 
singular point : 

a 


Y—Y 


Nar + f(y). (31) 


If 60, /6y+ 0, a can be so chosen that f(y) has a finite value. The electric field shows 
the same variation with the coordinate. The relative variation of the electric field 
E across the distance y then becomes 


AE Ay 


(32) 
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putting dy equal to the mean inter-electron distance and assume (31) to be valid up 
to this point. . 
_In an interval in the y direction, which under the present assumption is of the size 
of two times the mean distance between the particles, we cannot use the concept of 
a homogeneous region. Since very few particles exist in this region, we can assume 
at we have vacuum there. The amplitude of the few particles that do exist in this 
gion is limited by non-linear action, so that their influence may be neglected. The . 
dielectric constant on the two boundaries of the vacuum region now has the same 
_ absolute value but opposite sign, which means that the electric field retains its magni- 
_ tude but changes sign across the vacuum gap. In the vacuum gap we have a homo- 
geneous field. In that field those few charged particles present have to be treated as 
ics In the present treatment we shall neglect their presence. 
- According to the picture we have got we must, when calculating S, put o(y) =0 
J when y, —d<y<y, +6. We cannot expect any greater accuracy from such a proce- 
_ dure, since at the boundary of the region where we have assumed a continuous distribu- 
tion of charges this approximation does not hold. Nevertheless we may expect a 
_ qualitatively valid result. 
We introduce the quantity # defined by 


cosh }=6 {0% (-2)"l". (33) 


_Tf we let 6 tend to zero we obtain the solution in the form 


{Biome 


U Qaa2 2aar 2H 


i (34) 
S= a car . 
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We find the value of w2 for which the homogeneous equation (12) has a solution 
when we put the denominator of S equal to zero. An examination shows that there 
exists no solution for any value of w? when 0 < (2«/d)/(a?/b?) <3. This means that 
‘if the absolute value of the ratio of the electron space charge density and the back- 
ground space charge density is below a certain value we find no resonance. 
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Fig. 3 Fig. 4 


Fig. 3. Resonance condition when w?<(eB/m)? for a parabolic spatial charge distribution with 
stationary background space charge. 


Fig. 4. Resonance condition when (eB/m)? < w? <(eB/m)? + b? for aparabolic spatial charge distribu- 
tion with stationary background space charge. 


When (2«/d) (a#/b?) > 1 we have resonance for a frequency satisfying w? < (eB/m)?. 
Fig. 3 shows that the quantity (2«/d) (a?/b?) giving resonance tends towards 1, when 
the applied frequency tends to the gyro frequency of the electrons. When we increase 
(2a/d) (a?/b?) the value of w? giving resonance decreases, passes zero and decreases 
to — co. When the electron space charge density exceeds the background space charge 
density by a certain amount, our configuration is thus unstable. The limiting value 
depends on the ratio b/(eB/m). When this ratio is near zero the limit corresponds to 
(4/3)/(«/d)a? ~ (eB/m)*, which means that the electron density for the limit of the 
stable region is given by the gyro frequency of the electrons and not by the background 
charge density. 

Fig. 4 refers to the case when there is a singularity somewhere in the plasma 
region. Then the resonance frequency satisfies the relation (eB/m)? <w? < (eB/m)? 
+ 6. Further 4 < (2«/d) (a2/62) <1. 

Dissipative medium. If we allow collisions between the particles, the amplitude is 
limited everywhere in space. Generally the collisions are taken into account by the 
introduction of a friction term in the equation of motion. If the losses are sufficiently 
large we can use the concept of a homogeneous medium throughout space. Let us 
now assume a frequency corresponding to case 2 as introduced in the beginning of 
paragraph 2.1d. If we then let the losses tend continuously towards zero we find—in 
the limit—another result than the one we have just presented. If, however, we remem- 
ber that the simpler way of solving the problem may be used only as long as the rela- 
tive variation of the electric field is below a certain value throughout the plasma 
region, we find that, when we decrease the losses of our medium, we eventually arrive 
at this limiting value. Following these lines of thought we obtain the solution corre- 
sponding to the one found in literature, when the losses are sufficiently large, and the 
one we have presented in this paper, when the losses tend to zero. We thus find the 


same solution if we first introduce losses and let them tend to zero as if we start with 
the assumption of no losses. 
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(e) A neutral and homogeneous medium with several species of particles 


_ If we restrict the discussion to a homogeneous medium containing several species 
of particles we will find that all r are independent of the space coordinate. This means 
that they can be put outside the integral in the right member of equation (12). If we 
_ also assume that the medium is neutral and confined to a layer of thickness 2«, equa- 
- tion (12) becomes 
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The sum > a, 27, 
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is the same in all equations. If we transfer the sum to the right member, this 
member will thus be identical for all equations. Let us denote it by_S. The 


- solution is then 
: oat si 
teks { (2) -o| ; (36) 


My 


- Introducing these values in the equation defining S we find 


, U a -1 
ay : (22)'—o : (37) 


In Fig. 5 the denominator of S is given as a function of w*. We there have an asymp- 
tote at the gyro frequency of each species. We can also see that the number of zeros 
is equal to the number of species. But this is equal to the number of roots of the equa- 
tion. All roots in w? are positive and hence we can conclude that all values of w that 
make S-! = 0 are real in the case of a homogeneous neutral medium with any number 
of species of charged particles and a homogeneous magnetic field. At these zeros, S 


Fig. 5. S—1(?) for a neutral and homogeneous 
medium with three species of particles. 
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Fig. 6. The equivalent circuit for a neutral and homogeneous me- 
dium. 


becomes infinite and consequently all the r, become infinite. They therefore represent — 

resonances of our system. When U = 0 we can have a finite amplitude of oscillation 

at these frequencies. These oscillations thus represent natural modes of our system. 
At the gyro frequency of the wth species we easily find 


U —2 
alle BA ad sthnggeatlay 22 (38) 
m 
t=O when v+uy. F 


We thus see that when w is equal to the gyro frequency of the uth species the ampli- 
tude of motion for this species is finite and non-zero but zero for all others. 

Since we have found that all the natural modes of our system are oscillatory we 
must conclude that the equilibrium of our system is stable for perturbations of the 
kind we have discussed. In the special case when we let the magnetic field strength 
tend to zero we find 


ta OG (x a—o#). (39) 


We find that we have only one resonance frequency and that that is a real quantity. 
The natural mode of oscillation is thus also in this case oscillatory and our system is 
stable. The resonance frequency is the root of the sum of the weighted squares of 
the so-called plasma frequencies for the different kinds of particle present. The weight 
depends on the geometry. 

Our discussion shows that the introduction of a magnetic field splits the single reso- 
nance which exists when the magnetic field is zero into as many resonances as there 
are kinds of particles present. All these frequencies, except the highest one, approach 
zero if we decrease the magnetic field to zero. 


The admittance. The admittance is 


“i 2a/d 2, B\? e 
Tail Sie <<) with o=>a? (22) -o| : (40) 


Some simple calculations give us the circuit shown in Fig. 6. The contribution of each 
species is represented as a parallel resonance circuit with the resonance frequency equal 
to its gyro frequency. All circuits are coupled in series and in series with them all we 
have the capacity C, which, however, is negative. Because of the variation with fre- 
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From these expressions we find, in the same way as in the previous example, 


4 r = { (%2)"— 20 _ ae)" 
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U a; aS (42) 
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_ If we introduce {(e,B/m,)? — e,Q/m,€,}* instead of the gyro frequency in the previous 


/ case we can even now use the figure corresponding to the neutral case. We then find 
- that the asymptotes for those species that have a charge opposite to the excess charge 
_ always exist. The corresponding frequencies increase when we increase the absolute 
value of the excess charge. The frequency of the asymptotes for a species having the 
same sign of charge as the excess charge, decreases with increasing absolute value 
of the excess charge. This means that, for a certain value of the excess charge, the 
_ asymptote reaches the frequency zero. For a larger absolute value of the excess charge 
the corresponding asymptote does not exist any more. Each time this happens 
_ when we increase Q, the equation S—! = 0 loses two real roots and two imaginary roots 
_ appear. 
The disappearance of the asymptote occurs first for that particle with charge of 
appropriate sign which has the smallest value of e/m. 

Hence we have found that our system is stable if the excess charge does not ex- 
‘ceed a certain limit (but that it is unstable if the excess charge exceeds this limit). 
The limit is different for different sign of the excess charge. We also find that the 

range of stability cannot be increased by introducing new kinds of particles or by 
changing the densities of those already present assuming constant excess charge. Of 
course, we can decrease the range by introducing new kinds of particles. 


(g) The general case 


We shall study the general case of an inhomogeneous plasma consisting of several 
_ species of particles. We also allow it to be non-neutral. The last term and the two sums 
are the same for all the equations (12). Let us call the sum of all these terms S. Since 
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In equation (44) the two last rs independent of y. Let us call their sum 


k, We then find a com 4 
S=k {1+o}. — (45) 
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If we introduce the expression (45) for S in (46) remembering that & is inde- 
pendent of the space coordinate we get the following expression for & 


ee gee ee ee 4“ 
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Hence our solution can be written as 


1 (eB See |" at 2% iiss) (-4) 
=| (=) mene |) Rhea laa) Nie 
Thus r is expressed as a product of three factors. The first one is specific, whereas 
the second is the same for all species and gives the common space dependence, The — 
third factor is also the same for all particles. It is independent of the space coordinate. — 
When the first factor becomes infinite for some species, the second vanishes in 
such a way that the product is finite and non-zero. This means that the amplitude — 
is finite and non-zero for this species but zero for all others. For a neutral plasma 
this happens at the gyro frequency and applies to the whole plasma region. When 
the excess charge density is independent of the space coordinate this phenomenon 
concerns the whole plasma region but appears at frequencies which are modified by 
the excess charge density. When the excess space charge density is a function of 
the space coordinate, this phenomenon occurs only at discrete value of the space 
coordinate. 
We shall discuss the third factor for the special case of a neutral plasma, assuming | 
the following space dependence ~ 


2 
eu Y 


My Eo 
We then find that the denominator, D, of the third factor can be written 


2, ; 2 
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(a) a of eanennhs (51), (52) and (53). 

ae zero of D as a function of 2a/d. 
‘parameter u of equation (50) as a function of frequency o. 

@ D as a function of @ in the case when ads are resonances. In the dashed portion = 

is a singular layer in the plasma. 

‘ _ (e) Dasa function of in the case when ae are no resonances. In the dashed portion there 

z is a singular layer in the plasma. : 
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< ; in the earlier discussion, we have three cases 
_ (1) u>0. We introduce sh o=(u)? and find 


2a 2a 26 ¥ 
- ee i ida hes el 
(2) —1<w<0. We introduce sin y=(—w)? and find 
; OL Be) 
4 2 2) 


(3) w< —1..We introduce cosh #=(—w)?. If we now, as before, exclude a 
_ neighbourhood of the point where the integral becomes infinite, we find after 
; _ integration 


a O ot ioe 20 
ee = — ee 
Bris d d sinh 20° 58) 


; In Fig. 7 a we have given D as a function of uw. We find that D has one zero when 
2a/d>4 but none when 2¢/d<}. In Fig. 76 the zero u, is given as a function of 
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For the case of two species, Fig. 7c gives u as a function of w?. The asymptotes 
are found at the gyro frequencies. Using Figs. 7a and 7c we can derive DY as a 
function of w?. Fig. 7d illustrates the case when there are resonances, and we find that 
there are as many resonances as there are kinds of particles. Fig. 7e illustrates the 
case when there are no resonances. The resonances have been replaced by peaks of 
finite amplitude. 


In the last two figures there are points of discontinuity. They occur at the fre- 


quencies that give infinite amplitude in the densest part of the plasma region. If we 
decrease the frequency from a point of discontinuity, the region where we have 
infinite amplitude moves towards parts where the plasma density is smaller and, at 


the gyro frequency, reaches the edge of the plasma region. Consequently we find © 


resonance of a maximum value of D-! in the frequency range giving infinite amplitude 
somewhere in the plasma region. Resonance cannot occur when the amplitude is 
finite throughout the plasma region. 


(h) Discontinuous density distribution 


It is generally quite complicated to discuss the electromagnetic behaviour of an 
inhomogeneous plasma region. Therefore one is inclined to approximate the actual 
density distribution by a number of steps, each with constant density. In order to 
investigate the consequences of such a procedure we shall here discuss some calcula- 
tions of this kind. For our purpose we may treat the case where there are only two 
species, the same in all regions, assuming neutrality in each region separately and 
also that there is no magnetic field. The number of homogeneous regions is arbitrary. 

If we—in equation (12)—exclude the term related to the magnetic field and the 
one related to the excess charge density we find that the second derivative of the 
position vector has the same value for all species. If we assume harmonic variation 
in time, also r must have the same value for all kinds of particles at the same place. 
This also means that we can use one and the same equation for them all. In our case 
the number of equations is reduced by half. We get one equation for each homogeneous 
region. The right member of equation (12) has the same value in all these equations. 
Let us denote it by S. Letting the index of r refer to the region where the particle is 
situated, we obtain from equation (12) 


f, + O2r,=8, (54) 
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Here er and my, refer to one and e; and m,, to the other. If the phenomenon 
is harmonic with the angular frequency w equation (54) gives 


Ty = S (Q?2 — w?)"}. (57) 
Introducing this equation in (55) we find 
erage ob) Pa A 
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Fig. 8. The frequency characteristic for_a plasma con- 


In Fig. 8 the value of S is given for the case of three homogeneous regions. We 


_ easily find that there are as many asymptotes as homogeneous regions or equivalently 


as many resonances as regions. We also see that S vanishes at each of the frequencies 
defined by (56) but the denominator of (57) becomes zero for one of the regions. While 
the particles in the wth region have a finite amplitude at the frequency Q,,, the ampli- 
tude vanishes in the other regions. 

Thus, if we try to approximate a region where the space derivative of the charge 
density is continuous and non-zero everywhere in this way the number of resonances 
equals the number of steps introduced. Since the previous discussion has shown 
that such a system has only one resonance we must reject this “approximation”. 


2.2 Introduction to the case of finite temperature 


(a) Only one velocity represented in the unperturted case. Neutralizing 


_ stationary background 


A non-zero temperature is generally introduced by means of a pressure term in 
the macroscopic equation of motion. In the present paper we shall instead generalize 
the results of the previous chapter by introducing an unperturbed motion of those 
particles constituting the plasma. The development is not carried beyond the initial 
stage where the velocity spectrum consists of one discrete line only. 


We have shown that the motion of a charged region of one species is independent 
of the charge density distribution if the distance between the two plates is infinite 
but the charged region has a finite width. When there are two particle species with 
different values of e/m the situation is different. Now the constitution of the plasma 
is of interest also in the case of infinite distance between the plates. The important 
difference between these two cases is that in the first the particle velocity is a single- 
valued function of the space coordinate but in the latter the particle velocity has two 
values at every’ point in space, one for each species. The important factor is now the 
differential velocity. 

In the unperturbed state we have only one kind of particle and that moves back 
and forth periodically in the direction in which a perturbing electric field is applied. 
Also now there are two values of the velocity represented at each point in space, 
which is true also if we disregard the unperturbed part of the velocity vector and 
look only at the perturbation in the velocity. This means that the discussion of a 
charged region with finite temperature belongs to the second group of phenomena of 
the previous discussion. 
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As before, we shall introduce only one new feature at a time. We first assume that 
we have neutrality when there is no perturbing field. The neutralizing field may be 
due to a stationary space charge. Then, the unperturbed motion is linear in time inside 
the charged region. We shall simply assume that the particles are reflected at y =0 


and y=d. We also assume that the unperturbed absolute value of the velocity is — 


the same for all particles, and an equal probability for all positions between y =o 
and y =d. Finally, we assume that the perturbing electric field in the y-direction is 
obtained from plates at such great distance from each other that the distance can 
be treated as infinite. 


Integration of the equation of motion for an individual particle would be very much ~ 


simplified if we need not consider the reflection of the particles. This can be achieved 


if we introduce a mathematical space which coincides with the physical one when ~ 


0 <y<d. The mathematical space from 0 to -d is obtained by reflection of the physi- 
cal space in the surface y = 0. By repeated reflections in the surfaces y = 0 and y =d 
we can correlate corresponding points in the physical and the mathematical space. 

In the region 0 to d the electric field has the same value and direction in the mathe- 
matical and the physical space. In the region -d to 0 in the mathematical space mo- 
tion in the positive direction corresponds to motion in the negative direction in the 
physical space. This means that the electric field in this part must have a direction 
opposite to that in the physical space. Thus the electric field in the mathematical 
space must change sign for each distance d passed. 

We have just found that the electric field must be an odd function of the y-coor- 
dinate and at the same time periodic in space with the period 2d. Assuming harmonic 
time variation, the most general field satisfying these conditions is 


F=> A, sin nase, (1) 
n=1 


where we have used the symbol F for the field, since we want to use E for the part 
of the field due to the potential difference between the plates. 
The equation of motion in the mathematical space is now 


my=eF, (2) 


where F is the value of the electric field at the perturbed position of the charged 
particle. In equation (2) m is the mass of the particle and e its charge. 

In an experiment we always have a continuous velocity distribution. To simplify 
the treatment we shall, however, assume that all the particles have the same absolute 
value of the velocity in the unperturbed state. Half the number of particles move 
upwards and half downwards. The density in the physical Space may be JN. 

In the mathematical space all particles move upwards. All particles have thus, in 
the unperturbed case, the same direction and absolute value of the velocity. The den- 
sity is } NV. If the particles move along the y axis, the position y, of a certain particle 
in the unperturbed case is given by 


Y=vita, (3) 


where v is the velocity of the particle, t the time and a a phase constant. This phase 
constant is necessary because at any instant all positions in space are equally probable. 
In the perturbed case the particles are situated a distance 7 out of the unperturbed 
positions. If y is the perturbed position we have 
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4 Y=Yotn (Yost). . (4) 
put we now introduce (4) in (2) we find 
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We assume nz7j/d< 4x. We would like to approximate the cosine of this quantity 
_by 1 and the sine by the argument itselt. This is possible only if we restrict n to a 
finite region. In this section we shall formally introduce an unlimited series in n as 
_we know that in an experiment a very limited number of terms are of interest. Looking 
at the treatment from this point of view our procedure can be justified, even if it is 
_not fully correct from a formal point of view. With this approximation equation (5) 


becomes 
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This is an inhomogeneous linear differential equation with periodic coefficients. If 
we keep A, sufficiently small, ie. we treat only small perturbations, we can, as a 
first approximation, put the coefficient of 7 equal to zero and obtain an equation 
with constant coefficients. If we then introduce (3) into (6) we find 
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‘If we integrate and then introduce y, again we find 
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Here we would again like to neglect the terms na7/d. The same difficulty as above 
appears here and we treat it in the same way. We find 
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There are two sources of the electric field. The first is the potential difference bet- 
ween the plates, which are at infinite distance in the physical space. The other source 
is the space charge distribution. go sagl 

The first source creates a homogeneous electric field # = H, e o*'in the physical 
space. In the mathematical space it is a square wave with the period 2d. Hence it can 
be written 


= 4 
p=" H,é@? > (2n+1)~ sin (2n4+1) . (10) 
IU n=0 
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space. Furth ee perturbation for those particles moving upwards ir 

space is the same in the physical and the mathematical space but the pertur 
for those particles which move downwards in the physical space has c ifferent s 
in the physical and the mathematical spaces. If we then also remember that the d 
sity in the physical space is NV but in the mathematical space N we easily find the 
space charge contribution to the electric field at the position y . 


(4 N ele)n(—y)-(BNele)nty). (11) 
Thus the total electric field is 
F (y) =H + ($N e/eo) 4 (—y) — (3 N e/eg) n (y). ~ (12) 


If we introduce the expression (10) for H and (9) for 7 in the relation (12) we can 
compare this relation for the electric field with the expression for the field (1) introduced 
in the beginning of this section, The comparison between equation (1) and (12) gives 
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The unknown coefficients A, can be found very easily as each equation contains only 
one unknown. The equations for the coefficients with even indices do not contain 
the applied field # and therefore correspond to modes which cannot be excited by 
fields of the type discussed. The coefficients with odd indices depend on the applied 
field and can thus be excited. 

If we solve equation (13) for A,,,,, we find 


4 #F mV ~ Vv pet 
Aas 2 98g [ta [en+ + of -108} n+ 2-9 | (16) 


This is a particular solution for the equation when £ is different from zero, provided 


that none of the denominators vanish. If one of them does vanish we cannot find 
a stationary solution. 


The normal modes. Even in the case where we do not have an applied field #, equa- 
tions (13) and (14) may have solutions, but only if the following relation is satisfied 


a pat "i 
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Fig. 9. Graphical representation of equation (19), especially useful when the frequency of oscilla- 
tion is given. 


Fig. 10. Graphical representation of equation (19), especially useful when the plasma density 
is given. 


Let us introduce the notations 


m2 y2 2 2 
r= 302 /7 and s = 2 oe (18) 
Equation (17) can then be rewritten as 
(so )= 
2 
* Sn a 


This is the equation of a hyperbola. Only the branch for which r>0 is of 
interest here. It is reproduced in Fig. 9. 
The figure shows that when 


Gye! (20) 


equation (19) allows solutions for s <0, which means that we have natural modes 
with exponentially increasing and decreasing amplitude. Thus, if relation (20) is 
satisfied for integral values of n different from zero, the plasma is unstable. The num- 
ber of unstable modes equals the nearest integral number below Qd/zv. In the present 
case we have assumed that the plasma is confined by reflecting walls to the volume 
between the planes y =0 and y =d, and therefore the instability cannot cause an 
increase in volume. Thus the instability can only cause a change in the velocity distri- 
bution. We find that an one-velocity distribution is unstable if the common velocity 
in the unperturbed case is such that the time of flight of particles from one wall to 
the other and back again is larger than the period 27/Q. Our theory is too simple to 
tell what velocity distribution is necessary to give stability. However, we may state 
that, if the velocity of the particles is so large that the relation (20) is not satisfied 
for any integral value of n, the plasma is stable for the kind of perturbations discussed 


here. 
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Let us assume that w, v and d are given and that we want to know the plasma den- 
sities corresponding to natural modes. In Fig. 9 we can see that the density is a single- 
valued function of the parameter n. If we first discuss periodic disturbances we can 
distinguish two cases. The first is characterized by w >z0/d, i.e. s > 1. Fig. 9 shows 
that we then have a limited number of modes for which Q decreases with n?. Apart 
from these, there is an infinite series of modes for which the plasma density is a mono- — 
tonically increasing function of n?. ; 

‘The second case for periodic disturbances is characterized by w <zv/d, ie. s< 1. 
Fig. 9 shows that the plasma density corresponding to natural modes is then an infi- | 
nite series, the terms increasing monotonically with n?. 

We now come to the case where we have exponentially increasing or decreasing © 
perturbations, i.e. s<0. For rapidly growing solutions, i.e. for s < — n?, we haveno ~ 
natural modes for any value of the charge density. For solutions growing slowly 
enough we have an infinite number of solutions for integers n satisfying n > (—s)?. 

Let us now assume that ©, v and d are given. Then, for each value of n, equation 
(19) gives two values of w, one value corresponding to the branch s/n? <1 in Fig. 9 
and one value corresponding to the other branch of the curve in the same figure. 
This is seen more clearly in Fig. 10. We find that both values of w increase with n 
as long as we consider only oscillating solutions. 

In the domain of the complex plane for w obtained when excluding the axes, we 
find no solutions corresponding to natural modes. 

The impedance. Let us now assume that we have two plane metal plates at the 
finite distance D and that the region containing the charges is inside the smaller 
region d (d < D). In the region where no charges are present we then have the field 
EH, which we introduced previously. In the charged region we have the field F. 
The potential difference between the plates is given by E- D with any specified rela- 
tive accuracy provided D is sufficiently great. In the present case when D has a 
finite value the distribution of the charges also contributes to the potential difference. 

The potential difference across the plasma region U, is 


d 
U,= -| Fdy. (21) 


0 


But we have found that (see equations (14) and (16)) 


4 
Aon=0 and Aanni => 55 (1+ 8) (22) 
with 
2 2]: + 
o=02 for fensnF] |([ofenen te} ] -axfors {ens ®t ]) - 


(23) 


Introducing this expression in the definition of U,, remembering that F' is defined 
by equation (1), we find that the first term in A,,,, gives the potential difference 


— Ede’, which is due to the applied field #. The other term is due to the space 
charge. 
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: After performing the integration we find 


ey ie Pee ee $ si tee (24) 
I? nao (2n+1)2 te 
aif we write the terms as a sum of partial fractions we easily arrive at the Riemann 
-_£-function and the ‘Y’ functions in performing the summation [23]. Instead of the ¥- 
“functions we can introduce trigonometric functions. Finally, we arrive at 


O2 
A U,= —Bd-4Ed [H+ r— m)s (25) 


3s—a. tg ta12% 
(r2+47rs)* Laayst * 


4 where %2=r+st(r?+4rs)* and pis=} (26) 


_ Observing that, in the physical space, the field outside the plasma region is E we 
find, when the distance between the plates is D, the impedance of the system 


ge: ipl Sm, 5 tena 
wet i sage reps (=) aw? — (2 aie DE (27) 
d 


Here, the factor w? —Q2 appears in the denominator and it is of interest to determine 
_ the value of Z when w? =Q?. We easily find 


; ends oe 
tQOe, 202, 9 


2 a + 
pode 498 4 (67) 


Z 
+ 4a (6r)? 


when w?—Q?=0. (28) 


Consequently the impedance has a finite value at w? = Q2? if (6r)? =2y +1 where v 
is an integer. 

Because of the factor tg $22? in the expression for the impedance, the impedance 
becomes infinite when ~ 


xt=2y+1, (29) 


where » is an integer. To obtain the resonance frequencies for a given plasma density, 
we have to solve equation (26) for s and then introduce equation (29). We find the 
resonance condition 


= T (72 2 
s=a2tr4 (r2+427r) 4 (30) 


x=(2y+1)?. 
Equation (30) gives the same solutions as equation (19) when 1 is odd. The values 
given by equation (19) when 1 is even, are not represented in the present case. Con- 
sequently, we can use Fig. 10 in the present case if we exchange n? for x. The lower 
branch (Fig. 11), reckoned from the vertical tangent of the curve, then represents 
the solution obtained using the minus sign in equation (30). The rest of the curve 
represents the solution of equation (30) with the plus sign. We can conclude that the 
frequencies of the odd natural modes appear as resonances in the impedance of the 


system. 
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Fig. 11. Graphical representation of equation (30). 


As an illustration we give, in Fig. 12, the resonances for the special case when 
(2r)! =14. On the central scale we give the angular frequency divided by Q and on 
the lower scale the angular frequency divided by zv/d. With these scales as reference, 
the upper line shows the positions of the resonance frequencies. We have plotted 
one series of resonances above the line and one below. The first series is restricted 
to the range w > Q. The other begins with seven resonances for imaginary values 
of w. Then come five resonances for 0 < w <Q. Excluding the first thirteen resonan- 
ces, all the rest have w >Q. 

In any experiment there are losses. This means that the resonances are damped 
and the damping is more severe the higher the order of the resonance. In the present 
example it should be easiest to find the first seven resonances above the frequency Q, 
shown above the line in Fig. 12. The resonances shown below the line should be as 
strongly damped as the following resonances above the line—at least. The lower 
measuring frequency would surely make the damping even more severe. 

In an experiment it is usually more convenient to keep the measuring frequency 
constant and vary the charge density. The asymptotes of Z still correspond to 2,2 = 
(27 +1)2, 


_ [sear YP 
s+ (2+ 1)? 


Qr (31) 


Again we see the similarity with equation (19). The resonances coincide with the natu- 
ral modes of odd order. 


Equation (31) shows that the plasma density for resonance is a single-valued func- 
tion of the parameter y. 

If we order the resonances according to increasing values of », Fig. 13 shows that, 
to begin with, r decreases with y, i.e. the plasma density decreases for increasing ». 


a ; Ya, 5 - at 9 nN Lee ny ee 23° 255 ¥ 

15 17 19 21 eRe sy 29 31 33. «35 «637. C39 GI 43 .45.>N, 
-——____ 1 hes aes 1 dees 4 me — ee | wnilline 1 = ue oe os $$ | dl Ww 
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Fig. 12. The resonance frequencies when the resonance conditions are given. 


192 


ARKIV FOR FYSIK. Bd 19 nr 13 


= 
a 
ca 
i] 
= 
ot 
ne 
pid 
N 
w 
i) 
a 
nN 
x 
i) 
o 
w 
e 
w 
o 
w 
a 
w 
<] 
w 
o 
= 
< 


Vv 

ee ee Q 
1 2 o 
10 20 3 Ww 

d 


ge: 13. The resonance values of the plasma density when the frequency of perturbation is given. 


a 
_The number of resonances with this property equals the smallest integer larger than 
_4(s* — 1). Consequently, no resonances of this character exist if s <1. Next, the reso- 
-nances fall on the part of the curve in figure 9 with negative slope. Now the reso- 
nance plasma density increases with the order of the resonance. 
_ The total number of resonances for which Q <w equals the smallest integer larger 
than $ [(3s)* — 1]. No resonance with this property exists if 3s <1. 
_As an example we choose 


beg | oe 14 
s @ i 3 


{4(st-1)=63; 4((38)?-1)~11.614}. 


The resonance of lowest order appears for a value of Q slightly smaller than w. 

Seven resonances of consecutive orders appear for decreasing values of Q. The next 

“five appear at increasing values of Q but still have Q <w. Thus, altogether twelve 
_ resonances have the property Q <w. 


3. Cylindrical geometry 


3.1 Introduction 


It is well known that a homogeneous plasma cylinder of circular cross section has a 
resonance that can be excited by an oscillating electric field perpendicular to its 
axis. At resonance the reflected wave is identical with that reflected by a metal cylin- 
der with a diameter of roughly 2/3 of the wavelength. This is true whatever the 
diameter of the plasma cylinder. Since the value of r(é.B/é7) B just inside the boun- 
dary determines the reflection properties, our statement is true also for an inhomo- 
geneous column provided the density is a function of the radius vector only. We then 
take for granted that the treatment of the singularity is such that we find no dissipation 
in the plasma region when the collision frequency is zero. 

For a column diameter much smaller than the wavelength it is possible to approxi- 
mate the solution outside the cylinder by the static field. We then obtain infinite 
amplitude at resonance. Despite this defect our solution will not differ much from a 
more stringent solution. In the following we shall use this simplification and turn 
our attention first to the cylinder with elliptic cross section and then to the problem 
of an inhomogeneous density distribution. Because the numerical calculations in 
the latter case are very time-consuming, we have only brought those calculations, 
which are presented in § 3.3, so far that we see that there is only one resonance. 
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Fig. 14. Elliptical eylinder-coordinates. 


~~ 


3.2. Scattering by a plasma cylinder of constant density and elliptic cross section 


We shall assume that we have a plasma cylinder of elliptic cross section with con- 
stant density inside and vacuum outside. This cylinder is placed in a plane electro- 
magnetic wave propagated in a direction perpendicular to its axis and with the electric | 
vector perpendicular to the axis. If the cross section were circular we should expand 
the plane wave field outside the boundary in multipole fields and then treat the scatte- : 
ring of each component separately. We proceed similarly in the present case. ; 

The complete solution of the problem requires much tedious work. In order to. 
survey the problem we shall assume that the major axis of the cross section is much 
smaller than the vacuum wavelength of the exciting field. Thus we change the radia- 
tion problem into an electrostatic one, which means that we have to solve the Laplace 
equation instead of the wave equation. Our problem can be stated in the following 
way. We have a homogeneous dielectric rod in the form of a straight cylinder with 
elliptic cross section surrounded by an isotropic medium. At the boundary we prescribe 
an exciting field of multipole character and wish to find a field that is continuous 
everywhere inside the cylinder and falls off outside at least as r~, where r is the distance 
to the axis. At the boundary the potential and the normal component of the electric 
displacement must be continuous. 

The mathematical treatment -will be simplest if we use elliptic coordinates, &— 
and 7, which are related to the ordinary cartesian coordinates x and y by the defining 
equations 


x=c cosh € cos 7 (E> 0), i 
y=c sinh & sin 7 (0<n<2z). a) 
As seen in Fig. 14, € corresponds to the radius vector in ordinary cylinder coordinates 
and 4 corresponds to the angular coordinate. We also see that the curves, € = const., 


correspond to confocal ellipses, tending to a line of finite and non-zero length when 
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tends to zero. By a suitable choice of c in equations (1) the boundary of the dielectric 
? cylinder is made to coincide with a curve é = const. 

Let us denote the scalar potential by ®. The electric field is then — grad ©. We 
shall assume that the dielectric constant is a tensor (e) so that the electric displacement 
_is (e)E = — (e) grad ®. 
Z Since div D=0, we find 


+ div {(e) grad O}=0. (2) 
We easily find that this can be written 

z é e® o® é a® a® 

‘te mae Wa) Gene 0 as in. a ae = 

: se (se Sa") etm (oe* on) ,n Be 


5, where subscripts € and "indicate the € and 71 components of the vector re- 
_ spectively. If there is a magnetic field in the ¢ direction but no losses we can 


write [22] 
7 a b 

@=(_5 gi (4) 
Then equation (3) can be written 

2O &O_ 

ae + one Os (5) 


Attempting a separation of the variables we put 
® =u (&) v(n). (6) 


Since uw is a function of ¢ alone and v is a function of 4 alone equation (5) 
is split into 


on 


u’—n2u=0, 
(7) 


vy’ +n2v=0. 


Here n? is a separation constant. 


Since our problem is a physical one we must require that the solution is single val- 
ued. The defining equations (1) show that when 7 has increased by 2 7 we arrive at 
the point where we started. We thus require that the solution v(7) shall retain its 
value when 7 increases by 2 z, i.e. n must be an integer, n = 0, 1, 2,3.... 

We require a solution which is continuous and has a continuous derivative across 
the line = 0. This condition is fulfilled if we require the functions u and v be both 
even or both odd. Hence we get two pairs of solutions 


®,=cosh n€cosnyn and @®,=sinh né& sin 7. (8) 


Outside the cylinder we require a scattered field that decreases with distance. Such 


fields are ; 
®,=e "Fcosnyn and @,=e "* sin ny. (9) 


The potential of the exciting fields can be expanded in terms of 


r” cos n(y— a), (10) 
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where r and 9 are ordinary cylinder coordinates and « is an arbitrary angle corre-_ 
sponding to a rotation of the field pattern. The expression (10) can be rewritten as 


{ 
r” cos n(p—a) =r" cos NM cos Na+?” sin ny sin Na. (11)_ 
For a homogeneous field, equations (1) and (11) give 
r cos (py—«)=c cos « cosh € cos y+ ¢ sin & sinh & sin 7. (12) 


This equation shows that the two elementary waves in elliptic coordinates vary sinu- 
soidally with the angle «. Since € is constant on the boundary the maximum ampli- 
tudes of the cos 7 and sin 7 terms are unequal. 

The change to elliptic coordinates for higher modes may be made in the following 
way 


r” cos ng +ir" sin ny =(re'?)"=[c cosh (E+¢n)]” = (c/2)" (e&*'" + e--'")” (18) 
in particular, for n=2 and n=3. 
r2 cos 2m = $c (1+ cosh 2€ cos 27), 
r sin 2p=43c? sinh 2€ sin 2n, 
r= cos 3y= $c? cosh 3€ cos 37 + 33 cosh € cos n, 
rs sin 3p=4¢ sinh 3 sin 3n + 33 sinh € sin 7. 


Thus, each simple mode in cylinder coordinates corresponds to a finite number of 
simple modes in elliptic coordinates. 

Let us assume that the half axes of the cross section of the dielectric rod are a and b 
(a > 6). If € =€, is the equation for the boundary we get the defining equation 


b 
a engh oa (14) 


The plasma cylinder is situated in a field of the type 
®=y cosh n& cos ny +6 sinh vé sin nN. (15) 


In the special case n = 1, the first term corresponds to an electric field with the field 
vector parallel to the major axis of the cross section of the plasma cylinder and the 
second term corresponds to a homogeneous electric field with the field vector parallel 
to the minor axis. Outside the plasma cylinder we get a scattered field of the character 


O=c,e "* cosny+s,e"* sin n7. (16) 

Inside the plasma cylinder the field is 
=a, cosh né cos nn +b, sinh né sin nn. (17) 
The potential © and the normal component of the electric displacement are 
continuous at the boundary of the plasma cylinder. If the (relative) dielectric 


constant of the medium outside the plasma is ¢ these two conditions give 
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Cn=[—{(a—e)v +B? cosh n&,} y+ Pa ouhae) cosh n&,} bd] Se”®, : 
- 8 =[ {w—(a—e) sinh n&, } by—{(a—e)u+0? sinh n&} 6) Se"* 5) 


with 

: S-*=(a tgh n&)+ €) (a coth n&,+¢) +2, 
4 u=(a tgh n&,+e) cosh né,, 

v=(a coth n&,+ «) sinh n&p. 


In our paper [22] we computed the elements of the tensor (e) obtaining 


cree (fw) 


ME | \ m 

2 2 el 

Bes (fio) 2* | (2) -o| 
™m ME |\m 


where we have to sum over all species. In deriving these equations we assumed that 
the charged particles were free and moved under the influence of the electrie field 


and a static magnetic field, neglecting the influence of perturbations in the magnetic 
‘field. 


(19) 


- (a) No magnetic field 


In the case of vanishing magnetic field the solution is simplified considerably. If 
the exciting electric field is an even function of the coordinates € and 7, the excited 
fields are also even functions and similarly for odd functions. In this case the dielectric 
constant of the plasma cylinder is a scalar constant, ¢;, and the solution is 


ges at sett 
On = .=— : 
" e cosh né, +e; sinh n&, e sinh n&, +, cosh n&, 
4 (e—e,) sinh 2né,e"* _ $(e—6;) sinh 2n€,e"® (20) 
"e cosh né, +e; sinh n&, ” e¢ sinh né, +. cosh né, 
If we assume 6 =0 and y=1. If we assume y =O and 6=1. 


In an experiment, the field outside the cylinder is most conveniently measured in 
ordinary cylindér coordinates. To compare our result with experiments we need to 
convert our solutions to these coordinates. 

Generally we have, if we introduce (a + iy)/c =u, in order to simplify the notation 


e7f- in — g~ arcosh 4 —y—(u2—1)}. (21) 


If we then expand the root in powers of u™* we find 


=e 4 = 5 ig’ 7-8 95 
etalk tint te +g tau t+... (22) 
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Taking the nth power of this expression we find for »=2, 3 and 4 
eM ty dat Swot eu t..., a 


eo BED = hy 8 + Bu? ate take Ot neste» 238) 
ACEI) SE 4 atte aC cua ee Bs tee raiser 4 
e AED — Py t+ buf + sw b+... 4 


We thus find that the scattered field contains an infinite number of elemen tary 
solutions in ordinary cylinder coordinates. 
For the exciting field 


r 
a cos (p—«), 


+ 


where we have introduced the major axis of the elliptical boundary in order to 
make the expression non-dimensional, the scattered field is 


e-e (at+b\* [bcosa _, bsin.& 23 
2 5) | sae oe istic dhs a) 


where e§*'” is obtained from (22) and the conjugate complex equation. 


In an experiment we generally keep the frequency of the exciting field constant 
and vary the density of the plasma. This means that ¢, is our variable and 
can assume all values from 1 to— 9. 

In Fig. 15 c, and s, are plotted as functions of ¢;/e,. We find that C, and s, vanish 
when the dielectric constant has the same value inside and outside the boundary. This 
is self-evident as the dielectric constant then has the same value in all space. When 
£; decreases c, increases and, for a certain negative value of the ratio ¢; /€,, C, becomes 
infinite. This indicates a resonance effect and there is only one resonance for each 
value of n. The other curve in Fig. 16 shows that the sin ny mode is analogous in this 
respect. The resonances correspond to 

Ej &j 
+ 3 ctgh n&é and 76 teh: 2 f,. (25) 


0 0 
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3 ee 


oe becomes infinite, the plasma Anne reflects like a metal cylinder 
ite conductivity. The corresponding values of c, and s, are 


’ Si en bo 4 1 F e2nso a=4 1 
(Calmeta= ps ndmetar= ——3—" (26) 
troducing these in equation (20) we obtain ’ 
a =~ SI ay/e5 a —1+¢8/e, 
ed - Cn (Cn)metai coth n E, +e; Pe $n = (Sn) metal t gh n E, +8; /& (27) 


In Fig. 16 we notice that the cos my mode is reflected more strongly than the 
sin 27 mode. From equation (14) we easily find 
(a +b)" —(a—b)" 


(a+ b)" + (a—6)” (28) 


tgh n 50 


q ‘When the ellipse does not differ much from a circle, (a—6) /(a+ 6) is much less 
_ than one and we find, approximately 


4 r b 

i tgh n€,~1-—2 (=). 

4 : (29) 
: a 

4 ctgh nee +2 (5 a)’ . 

3 : 

. When }/a is much smaller than one, the resonance conditions are simply 


_ &/&* —na/b for the cosny mode and n™ order led (30) 
e,/&)¥ —nb/a for the sin my mode and n™ order resonance 
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of n. We have found 


From equation (20) it is evident that ai 
710 a 
Cn (iy Eo) = — Sn (Eq, &i)s 8n (Eir Eo) = — Cn (Eq; €i)- (31) 


Thus the geometric means of the resonance values of ¢;/e, for the cos nn mode 
and the sin »7 mode is equal to minus one, i.e. equal to the value for a circu- 
lar cross section. - 


(&/€9)oosnn : (&:/€9)sin nn == 1}. ~ (32) 


Fig. 18 illustrates the same conditions as Fig. 16 but for negative values of e/e, 
and with logarithmic axes, The abscissa has been normalized so that the resonance | 
occurs at minus one. The abscissa minus infinity corresponds to metal reflexion. To 
the right of the resonance in Fig. 18 the curve approaches the resonance level of 
the sin ny mode. It is reached at the plasma frequency. The sign is, however, opposite. 

The corresponding curve for the sin x7 mode is a mirror image of this curve relative 
to the ordinate axis, normalized to the sin ny resonance. Metal reflexion corresponds 


e=06 
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Fig. 19. Distribution of resonances for three values of the eccentricity of the plasma cross section. 


200 


ARKIV FOR FysiIk. Bd 19 nr 13 


Table 1. 
: i a Sa Sa ROE saeeeae eee Eee eee 
4 fom bass) OLS 0.6 0,28 
—SSSSSSSSSSSSSSSSSSSE 
ee 
4 (a) Resonance for cos n7 terms. €,/&9 =f. (n, b/a). 
¢ % 1 — 1.2500 — 1.6667 — 3.5714 
2 — 1.0250 — 1.1333 = 1.9257 
4 3 — 1.0028 — 1.0318 — 1.4330 
4 — 1.0000 — 1.0078 — 1.2225 
5 — 1.0020 — 1.1194 
(6) Resonance for the sin n7 terms. E/E =f, (n, b/a). 
1 0.8000 0.6000 0.2800 
2 0.9756 0.8824 0.5193 
3 0.9973 0.9692 0.6978 
4 0.9997 0.9922 0.8180 
5 0.99997 0.9981 0.8934 


to the left infinity, whereas at the right infinity, ie. at the plasma frequency, the 
reflexion is like the cos ny mode for a cylinder of infinite conductivity, again with 
opposite sign. 
Another consequence of the symmetry of equation (31) is that, the reflected energy 
is equal for the two modes when ¢;/e, = — 1. This is easily seen in Fig. 17. 
In Fig. 16 we can see that the cos 7 mode reflects more energy than a metal cylinder 
~up to a value of ¢;/e, slightly below zero. We also see that the sin 7 mode reflects 
more energy than a metal cylinder even for a slightly positive value of the same ratio. 
_ For a circular cross section all orders of resonance occur at the same value, ¢;/e, 
= — 1. For an elliptic cross section each resonance splits into two, one above and one 
below e/e, = — 1. Further resonances of different order occur, at different values of 
éi/e,. Some examples are given in Fig. 19. Thus for an eccentricity of the ellipse of 
0.6, we distinguish the dipole, the quadrupole and octopole lines. The higher modes 
lie too close together to be shown in the figure. When the eccentricity has increased 
to 0.8 we can distinguish 8 lines and when it has increased to 0.96 at least 10 lines. 


Fig. 20 Fig. 21 
Fig. 20. Electric field lines of the cos 7 mode. 


Fig. 21. Electric field lines for the sin 7 mode. 


20 and 21 illustrate the field pattern for the cos 7 mode and sin 7 mode. For 
the lowest mode the field is homogeneous inside the plasma cylinder. When the eccen- 
tricity approaches unity the sin 7 mode obviously tends towards the plasma resonance 
of a plane layer. Equation (30) shows that this limit resonance occurs at ¢/e, = 0. 


(b) Non-vanishing magnetic field 
We shall assume that the applied frequency is much higher than the gyro 


frequencies for any ions present, so that the ions may be considered immovable. 
Introducing g=eB/m and s=tgh n&, we can rewrite equation (18) 
2 Q? : 2 
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The resonance frequencies are given by 


1 é e Pye 1 s 
nan5 (242048) + fo (2-2) Pa (s +e) (;+e)} (24) | 


and the zeros by 


1 2 2 
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We notice that the scattered wave is elliptically polarized, i.e. there is a non-zero 
field along both axes of the cross section with a phase difference of $7 radians. The 
two components generally have different amplitudes. 

The resonance conditions are illustrated in Fig. 22. As the surrounding medium is 
generally vacuum (air), we assume ¢ = 1. The figure shows that there are two reso- 
nances when q@ is smaller than the electron gyro frequency, but only when the gyro 
frequency exceeds w. In the well-known case of a circular plasma cylinder, the 
two resonances degenerate into one, which appears when the ordinate is 2. The devival 
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2 the point when the abscissa is zero. For a tube of greater eccentricity 
e lower resonance appears at a lower value of the ordinate which must however 
eed unity for no magnetic field; in the latter case the curve in the figure is a straight 
e. The lower resonance is thus limited to the small region between the parabola 


for increasing eccentricity and also for increasing magnetic field. 
_ The higher of the two resonance frequencies is always higher than that correspond- 
ing to the parabola of a circular cylinder. In Fig. 23 we see that the ordinate increases 
for increasing eccentricity and increases for increasing magnetic field. There is no 
upper limit. 
At this point it may be appropriate to stress that a real magnetic field has two ef- 
_ fects: one is the one we have explained above but the other one is the change in the 
distribution of charge. We must therefore be very careful in applying the above cal- 
culations to experimental data. : 


Fig. 23 Fig. 24 


3 Fig. 23. The amplitude c, (along the axis perpendicular to the paper) as a function of 02?/@? and 
i F g?/? for an exciting field parallel to the major axis. 


Fig. 24. The amplitude s, (along the axis perpendicular to the paper) asa function of Q2/w? and 
g?/@? for an exciting field parallel to the minor axis. 


203 


E. AsTROM, Electromagnetic behaviour of a plasma 


Let us assume a linearly polarized electric field parallel to the minor axis of the 
cross section of the plasma cylinder, which gives the resonances shown in Fig. 24. 
The scattered field along the minor axis has a zero for a density that is independent 
of the strength of the magnetic field. This zero coincides with the highest of the 
resonance frequencies when the magnetic field vanishes. Consequently there is only . 
one resonance when the magnetic field vanishes but otherwise two resonances with ; 
a zero in between. As the component of the scattered field along the major axis has — 
no zero, it has always two resonances. However, its intensity is zero for a vanishing 
magnetic field. i 

For an exciting electric field along the major axis there is a zero in the scattered 
field in the same direction and this zero is independent of the strength of the magnetic : 
field. Its value corresponds to the lower resonance without magnetic field. As before, — 
there is then only one resonance for vanishing magnetic field. For finite value of the 
magnetic field strength there are two resonances with a zero in between. There is no 
zero for the field component along the minor axis. However, this component vanishes 
when the magnetic field strength vanishes. 


3.3 The inhomogeneous plasma cylinder 


In all experiments on the scattering from a plasma cylinder, the plasma has been 
produced by a gas discharge. It is impossible to obtain a uniform radial distribution 
of the density in a gas discharge. Therefore, one would like to have a theory for the 
scattering from an inhomogeneous plasma cylinder. The main difficulty has been and 
is still the singular surface where the dielectric constant vanishes. To introduce a 
step function for the density or a collision frequency that is finally put equal to zero 
only serves to hide the difficulty. The first method leads to one resonance frequency 
for each step introduced. Actually there should be only one resonance. The second 
method leads to losses also in the limit when the collision frequency is zero. 

In order to resolve the difficulty, we must first recognize that the concept of a di- 
electric constant is not applicable at and near the singular surface. According to the 
simple theory, the electric field strength approaches infinity at the singularity. When 
defining a dielectric constant one assumes that the relative change in the macroscopic 
fields is small within distances typical for the microscopic structure of the medium. 
This condition is not fulfilled near the singularity and we cannot use the concept of 
a dielectric constant in this region. Thus we must ask: How shall we treat the problem 
at and near the singularity? One possibility is to introduce particles distributed at 
random in space with a certain mean density and take account of the fields created 
by each of them. This, however, seems to be too difficult to be feasible. 

For a continuous medium we find that A (0#/0r)/E increases to an infinite value 
at the singular surface (r is the coordinate along the density gradient and A the inter- 
electronic distance). When the change in the field in the distance A is comparable 
to the field itself, ie. when A (0H /0r)/H is of the order unity, the idealization of a 
continuous medium must fail. The region where the concept of a dielectric constant 
is invalid has a width of roughly 2A and thus contains very few particles. We may 
consider this region to be a vacuum as this cannot introduce any significant error. 
We can show that the position of the boundary between the plasma and the vacuum 
region is not critical. Despite the fact that such a treatment is far from strict, it 
will help us to survey the problem of the inhomogeneous plasma. 


204 


ARKIV FOR FYSIK. Bd 19 nr 13 


_ We shall idealize the physical problem and state it mathematically in the following 
way. Let us introduce two surfaces, one on each side of the singular surface, and at a 
distance A, from it. As ¢ is zero on the singular surface the dielectric constant has 
the same numerical value but opposite sign on these two surfaces. Between the sur- 
faces there is a vacuum and outside the plasma is undisturbed and assumed to be a 
_ continuous medium. The potential and the normal component of the electric displace- 

ment shall be continuous at the two surfaces. Then we let A, tend to zero. In the limit 
_ we obtain a well-defined mathematical problem. 


We shall assume a plasma in the form of a circular cylinder with the density being 

a function of the radius only. We also assume a plane electromagnetic wave in the 
‘surrounding medium, whose direction of propagation and electric vector are both 

perpendicular to the axis of the cylinder. This field can be expressed as asum of multi- 
pole fields. In this paper we shall investigate the influence on the plasma cylinder of 
the multipole fields. To simplify the problem we shall assume that the radius of the 
tube is much smaller than the wavelength of the exciting field. It is then possible to 
use an electrostatic approximation, thus: We have a dielectric rod of circular cross 
section and infinite length and with a dielectric constant which is a function only of 
the distance from the axis. It is placed in a multipole field. In this approximation all 
fields are assumed to be independent of time. Since Maxwell’s equations give 


curl E = 0, (1) 
we can introduce a scalar potential so that 
E = grad O. (2) 
| Further divD=0 and D=ez«,E. (3) 
From equations (2) and (3) we obtain 
div (e grad ®) = 0. (4) 


Tf we introduce cylinder coordinates, equation (4) is 


1 o aM\ 14 5) 
=0. 5 
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To separate the variables, put 
O = R(r) 0 (9). (6) 


Then equation (5) is split into equations (7) and (8), where n” is a separation 
constant 
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This equation has three singularities with the exponents — 
Se 
nm and —n atr=0, E 
0 and 0 atr=a/b, 
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p-52(n+3) at r=oco.| 


The solution can be expressed using hypergeometric series. At r =0 and at r =a/b 
the difference between the characteristic exponents is integral. Hence only one of 
the linearly independent solutions can be expressed as a power series. At these two 
points the other linearly independent solution contains a logarithmic term. At infinity, 
however, both the linearly independent solutions are expressible as power series 
mln. 

want a solution that is finite at r = 0 and choose the one without the logarithmic 
term. We then have a solution that is determined for r <a/b.. 

At r =a/b this solution increases logarithmically. For large values of r our solution 
can be expressed as a sum of two linearly independent hypergeometric series of the 
variable 1/r. Also this solution increases logarithmically when we approach r = a/b. 


As explained earlier, we now have to determine the coefficients of the two hypergeo- _ 


metric series in such a way that the two solutions, the one outside and the one inside 
the singularity, approach the same value when r tends to a/b. The derivatives shall 
approach the same numerical value but have opposite sign on opposite sides of the 
singularity. This can be written as follows 


lim (®,~a)044 — ®,-a-,a) =90, 
A->0 
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The singularity will have its standard position if we introduce the new independent 
variable 


oe = dria. (13) 


In order to get one of the characteristic exponents at the origin equal to zero, we 
also introduce a new dependent variable 


Pe Res, (14) 
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ce this solution in the hypergeometric series of the variable 1/(1— 0). 
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j In order to get a solution on the other side of the singularity we retain the 
ee me terms of the solution R=" F (a,b; a+, @) near the singularity. We 
then obtain 
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Near infinity the solution can be expressed in two linearly independent hyper- 
geometric series. This fact can be used to obtain a solution when g exceeds 


unity. We have 
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As mentioned before, p and q are to be determined from the conditions that 
R shall be continuous at @=1 and that the derivative of R shall approach the 
same absolute value but with opposite sign when we approach g=1 from the 
two sides. Thus 


lim {R(1-A)—R(1+A)} =0, 
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> Le. red /dr, is also continuous. If we then introduce the new variable 


(37) 


- know that it must be continuous at the boundary of the plasma cylinder. 
Using (9), (18) and (37), equation (10) can be transformed into a Riccati equation 
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We want a solution of (10) that is finite near the origin. If we write this solu- 
‘tion as a power series, retaining only the first two terms, we obtain 
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The boundary value problem Vole 4 Bian 


Inside a circle of radius ry the medium is | aéseribied by equation (9) and « ) 
side 7) there is vacuum. In the latter region we assume a source with /9 
and a field of the character R=«a(r)/r)", where « is an arbitrary constant. 
the field in the vacuum region is ine 


n n 
R-(*) +0(%)' P= Fe: 
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At the boundary the normal component of the electric displacement, ie. e(dR/d ”, 
or equivalently 7, must be continuous. The tangential component of the electric 
field must also be continuous, i.e. R must be continuous. : 

_ Evo L—o (t9/r)"” 


mee 1+ (r,/r)* when r> To: (46) 
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At the boundary we find 
N Er, +0 1 BAe, 
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where 7;,-9 is a solution of equation (38). Hence 
-1 
x=(1 4 Er=0 Mre-0 ) (1 _ &r=0 ‘=0) (48) | 
Erj+0 = Eri0 2 ' 
In the special case where the absolute variation of the dielectric pe 
from axis to boundary is much smaller than unity we can use an approximate | 


value of 7 at the boundary. Using the approximation (43) and remembering that i 
é=a(o—1), we find 
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on 2s case « is still proportional to the difference between the dielectric con- 
stant outside the cylinder and its value at the origin. Because of the factor 
4n+ 1 in the denominator, « is relatively insensitive to é;~o. 


— Resonance a 


Assume a plasma cylinder situated in a medium of dielectric constant ¢,,,o. 
_ Equations (9) and (13) give for the value of o at the boundary 


0 =1—€&,-0/ér-0. (53) 
7 At resonance we further require that « shall be infinite. Then equation (48) gives 


> 


be | 


= 


aie &r,+0/Er=0- (54) 
- Eliminating ¢,-9 between equation (53) and (54) we find 
n= (Er,+0/Er,-0) (1 —@). (55) 


_ This is a straight line as seen in Fig. 25. The dielectric constant just inside the boun- 

, dary is constant for each line but ¢--9 varies along the line. We assume that the 

_ dielectric constant of the outer medium, ¢,,,9, is positive. Equation (54) then shows 
that 7 is negative when the dielectric constant at the origin is negative. 


14:2 211. 


E. AstROM, Electromagnetic behaviour of a plasma 


The straight line on the extreme right corresponds to a plasma density equal to 
zero at the boundary. For increasing density at the boundary the line becomes steeper 
and the vertical line corresponds to e =0 at the boundary. In the next quadrant € 
is negative at the boundary. The cusped curve in Fig. 25 corresponds to the solution 
of equation (38). : ; 

For ¢,_, > 0, i.e. 7 > 0 there is no intersection between any of these straight lines 
and the solution of equation (38). Thus, when the dielectric constant is positive at 
the origin there exists no resonance for a cold plasma with constant density gradiant 
in radial direction, assuming the density to decrease outwards. 

We have assumed vacuum outside the cylinder. For vanishing plasma density at 
the boundary the resonance then occurs when g = 1.40, i.e. when the dielectric con- 
stant at origin is — 2.5, This intersection occurs at the minimum on the curve of 
the solution of equation (38). This means that for increasing values of the plasma 
density at the boundary the crossing occurs at greater values of 7 and equation (54) 
shows that ¢,_) must decrease (both 7 and ¢,_») are negative quantities). In other words, 
to remain at resonance when we increase the plasma density at the boundary we 
must also increase the density at the origin. When the dielectric constant has reached 
the value zero at the boundary the plasma density at the origin must be infinite. 

Here we must add that a cylinder with infinite density does not show any resonance. 
The point corresponding to ¢ = 0 at the boundary is thus a singular point which we 
can exclude from further discussion. 

When further increasing the density at the boundary, we have to decrease the 
density at the origin in order to retain resonance and ¢,, »= —1 corresponds to a 
homogeneous plasma (the intersection is at @ = 0). Up to this point the intersections 
have been in the right half plane in the figure, corresponding to a plasma density 
increasing inwards. If we increase the density at the boundary further, the inter- 
sections fall in the left half plane, corresponding to a plasma that is densest at the 
boundary. However, the density at the origin must be larger than that corresponding 
to'é,, = 0. 

Thus, for a given plasma density at the boundary we find one and only one value 
of the density at the origin that gives resonance. 

Let us now assume that the ratio of the densities at the boundary and at the origin 
is constant. As a measure of the density. we use the plasma frequency Q. If Q, is 
the value at the origin and Q, the value at the boundary we put 


Q, ee 
—- =k=positive constant. (56) 
Q, 


Then equations (53) and (55) can be written as 
@ = (1 — B).Q§ (Q5 — w)“}, (57) 
= &,40(1—@(1—k)~*). (58) 


We again find a linear relationship between 7 and 9. All lines pass through the point 
4 =Né,,.9 and g =0. When the outer medium is vacuum, this point corresponds 
to a plasma with the density zero. The parts above this point correspond to ¢>0 
in the plasma. Since these parts do not intersect the curve corresponding to the solution 
of equation (38), they are not shown in Fig. 26. The points corresponding to Q, =m 
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in Fig. 26 and is identical with the corresponding curve in the previous figure. 

_ When £ is positive and smaller than a certain limit (about 0.44 for n = 1 and 0.34 
for n = 2) there are three intersections. This means that if we increase the density 
_ while keeping k constant, we find three resonances. If k is above that limit there is _ 

_ only one intersection and hence only one resonance. 
_ The vertical line £ =1, corresponds to a homogeneous plasma. The lines to the 
_ left correspond to a plasma that is densest at the boundary. — 


— Vanishing “reflected wave”’ 


The constant « vanishes when 


Yq ey, Er,+0/ €r=0- (59) 


If we assume ¢,,,) >1 and ¢,_) <1 the region —n &,,;9 <7 <0 does not correspond 
to any physical situation of the types considered here. Eliminating ¢,_) as before, 
we find (see Fig. 27) 


Se ke as, 
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Tay (Er,+0/€r,-0) (o> 2) (60) 


4 when the density is given on both sides of the boundary. Under the conditions given 
~ here no solutions can be found. Under less stringent conditions solutions are possible. 
_ Similarly, we find (see Fig. 28) 


n=Nér,+0 {0 (1—#) 1-1} (61) 
213 


or 
} = cc? oY wae os 


he vaaee tl as 


<e ©, , @Y Sgative rs 8 + otter 
\ cP suet hie he (uae chee) Se Ins AF ma 
ey cto svaln © sauve Rb eh adie Miot duet gin il el s 
TA oes #h ie manor ia pipes cigie must be itil | 

Neri: 0) oder edie (onder nnttetees efi aioe: 


bind 0d. Cadigaie, That 4s ruil PMY idethy 
‘ 


oF 


ah ee 
a 
oa 


Fig. 29. The “‘reflection coefficient” for the dipole mode as a function of the dielectric constant 
the origin when the plasma density vanishes at the boundary. 


when the ratio of the density at the boundary and at the origin is given. There is no 
value of k such that the reflected wave vanishes. Under less stringent conditions solu- 
tions are possible. i 

If we plot « as a function of ¢,_) when the density at the boundary is zero and there 
is vacuum outside, the result is that given in Fig. 29. We find « =0 when fap = 0 
as this corresponds to vacuum everywhere. Then « increases when « decreases. (The 
derivative at ¢ = 1 is finite; the vertical tangent in the figure is due to the non-linear 
scale.) At ¢ =0 there is a discontinuity in the slope but « remains finite, non-zero 
and continuous. Resonance occurs for ¢~ —2.5 and when « decreases further, « 
decreases and finally approaches unity. { 


4. Measurements on plasma resonators 
L. Tonks [1] seems to have been the first to publish experimental evidence for 
the existence of several resonances of dipole character. Two peaks are clearly seen | 


in his curves. 


214 


ARKIV FOR Fysik. Bd 19 nr 13 


4 Tonks used a gas discharge in mercury vapour in a glass tube of circular cross 
4 section. The discharge tube was placed between parallel condenser plates, which were 
connected to a Lecher line. He measured the reactive and the resistive parts of the 
impedance due to the presence of the plasma as a function of the plasma density. 
- The measurement was made for different values of the temperature of the condensed 
_ mercury (from 13°C to 45°C), for different values of the frequency (from 159 to 366 
- Me/s) and for different sizes of the condenser plates between which the discharge 
_ tube was placed. Nine measurements are tabulated in the paper. Tonks mentioned 
that there is also a third resonance in the cases when “‘the applied field was decidedly 

_ non-uniform’’. 
The next paper of interest in this connection was published 20 years later by D. 
_ Romell [4]. He used a gas discharge tube of the same radial extension and filled 
with mercury vapour. The curves refer to a case when the temperature of the con- 
_densed mercury is about 25°C. The measurement was made at a frequency of 109 c/s. 

Romell used a radiating antenna at 1.5 m distance from the tube and measured the 

reflected signal at about the same distance. The published curve shows three reso- 
 nances. 

The third paper of interest was published by A. Dattner [5]. He used a tube having 

a diameter 1/5 of that used by the other experimenters and at the same time a fre- 

quency 5 times that used by Romell. Further, he placed the discharge tube in a rec- 

tangular wave guide, operating in the 7'E,, mode, in such a way that the electric 
_ field was perpendicular to the axis of the tube. He measured the energy passing 

the cross section where the discharge tube was placed. At resonance, only part 

of the energy of the incident wave passed the section containing the discharge 
- tube. In the published curve at least six resonances can be seen. Also in this case 
the gas was mercury vapour but the temperature of the condensed mercury is not 
given. 

Extensive experimental data have already been published, but more are needed 
if we want detailed information about the influence of the different parameters and, 
before a quantitative comparison between theory and experiment is possible, the 
theory has to be refined. 

Our theory shows that the product of the frequency and the distance between the 
reflecting plates is an important parameter if we assume that the electron temperature 
is constant. In the experiments performed by Romell and Dattner the mean free 
path of the electrons, corresponding to the temperature of the condensed mercury 
is several times the diameter of the tube. This means that we should expect a large 
number of resonances. In the other measurements the mean free path is only some- 
what larger than the diameter of the tube. Therefore, we should expect fewer resonance 
peaks, which agrees with the measurements. 

In Fig. 30 the abscissa is the product of the frequency and the diameter of the dis- 
charge tube. The ordinate is the ratio of the density at the second maximum and that 
at the major maximum. The measured points are grouped along the upper line. 
Along the lower line we have plotted the ratio of the density at the tertiary maximum 
and that at the major maximum. The result seems to indicate that the positions of 
the minor maxima are determined by the value of fd. Despite the fact that the gas 
density, the linear dimension and the frequency used have been varied roughly a 
power of ten, the spread of the points in the figure is not very large. It therefore 
appears that the influence of the mean free path on the positions of the mmor maxima 
is a second-order effect. 


, 


215 


Fig. 30. A collection of experimental data on the relative position of resonances as a function 
of the quantity fd. 


Qualitatively the result is what we should expect from our theory, but at present 
we cannot make a quantitative comparison. 

Tonks, 1931 [1] published some measurements on a tube of elliptical cross section. 
His curves show two resonances for the major axis and only one for the minor axis. — 
The experiment thus seems to indicate that when we deform a circular cross section 
each dipole resonance is split into two. 

Dattner has made some measurements on a plasma resonator with elliptical cross 
section. The tube was placed in a rectangular wave guide in such a way that the elec- 
tric field was perpendicular to the axis of the tube. The measurements were made at 
a frequency of 5 x 10° c/s. From the published data we can clearly see that the resonance 
pattern is different for the case when the major axis of the cross section is along the 
electric field and when the minor axis ofthe cross section is along the field. We can 
easily see that the plasma density necessary for the main resonance is larger in the first 
of these cases. 

The results of some measurements made with Dattner’s equipment are shown in 
Fig. 31. The abscissa is proportional to the electron density in the tube and the ordinate 
is proportional to the transmitted energy in the wave guide. When the angle between 
the electric field and the major axis of the cross section is zero or 471 we find the curves 
that Dattner reported. For intermediate values of the angle we find a weighted mean 
of these two curves. This may be interpreted thus. There is one resonance pattern 
corresponding to each axis of the cross section and the excitation of these two patterns 
is proportional to the component of the electric field along the respective axes. 

We have discussed the case of elliptic cross section and homogeneous density with 
the tube situated in free space. In the measurements the density was not homogeneous 
and the tube was not situated in free space. Therefore, a quantitative comparison 
between theory and experiment is not fully appropriate. Theory and experiment 
agree in that we have one set of resonances for each of the axes of the cross section 
and in that a larger density is necessary for resonance when the major axis is parallel 
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_ Fig. 31. The measured value of the energy passing the cross section in a wave guide containing a 
_ plasma column of elliptical cross section, as a function of the angle between the electric field and 
, the major axis of the plasma cross section. 


_ 
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_ to the electric field. The major axis gives a more pronounced resonance than the minor 
_ axis—as the theory predicts. The order of magnitude is correct. 

_ If we wish to make a rough comparison between theory and experiment we may 
take the measurements made by Dattner and those made by Romell on tubes with 
_ circular cross section. We may call the radius R. The mean distance between two points 
on the circumference is 4R/z. For a Maxwellian electron-velocity distribution with 


_ the temperature 7’, we find the mean value of |v|~? 


[ol t=(b2m/kT) +. 


- 


j 
: 
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Hence the mean value of the time of flight from one point on the circumference to 
any other point on the circumference is 


If we then take 7 =410°°K, R =3x10-*mand f = 5x 10°c/s asasuitable approxima- 
tion in Dattner’s experiment, we find s = 55 and 4(s* — 1) 27. This last quantity — 
should equal the number of resonances below the density making the plasma fre- 
quency equal to the frequency of measurement. The extrapolated experimental value — 
is about 8. ; 

In Romell’s experiment /R is roughly the same as in Dattner’s experiment and 
the observed number of resonances is 3. The extrapolated value of the number of 
resonances is here roughly 6. 

Despite the fact that our theory is not developed for the cylindrical case nor 
accounts for a continuous velocity distribution, we obtain the same order of magni- 
tude for the number of roots as can be extrapolated from the experiments. From 
this fact we cannot draw the conclusion that our theory explains the experiments. 
However, our results can be taken as an encouragement for further research along 
the lines followed here. 


4 
| 
Let us use this as an estimate of d/v in our definition of s (equation 18 of § 2.2a). — 
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